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Abstract 

We present a method for generating exact interior solutions of Einstein’s equations in the case of 
static and axially symmetric perfect-fluid spacetimes. The method is based upon a transformation 
that involves the metric functions as well as the density and pressure of the seed solution. In the 
limiting vacuum case, it reduces to the Zipoy-Voorhees transformation that can be used to generate 
metrics with multipole moments. All the metric functions of the new solution can be calculated 
explicitly from the seed solution in a simple manner. The physical properties of the resulting new 
solutions are shown to be completely different from those of the seed solution. 
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I. INTRODUCTION 


General relativity is a theory of gravity and as such should be able to describe the held of 
all possible physical configurations in which gravity is involved. All the information about 
the gravitational held should be encoded in the metric tensor which must be an exact solution 
of Einstein’s equations. In this work, we will focus on the study of the gravitational held 
generated by compact objects like planets or stars, which can be considered as independent 
of time. In this case, the problem of describing the complete held generated by the source can 
be split into two correlated problems, namely, the interior and the exterior rotating helds. To 
handle the corresponding held equations, one usually assumes axial symmetry with respect 
to the rotation axis. It is well known that the exterior held of an arbitrarily rotating mass 
can be described by the Kerr spacetime [1]. As for the interior held, the situation is more 
complicated. In fact, a major actual problem of classical general relativity consists in finding 
a physically reasonable interior solution for the exterior Kerr metric, which could be applied 
to describe the interior held of realistic compact objects. One usually assumes that the 
matter inside the object can be described by a rotating perfect fluid. Since the discovery 
of the Kerr metric in 1963, many attempts have been made to find the corresponding exact 
perfect-huid interior solution. It seems that the rotation parameter for a perfect fluid leads to 
important difficulties in the context of Einstein’s equations. In this work, we will concentrate 
on the problem of perfect-huid solutions, without considering the rotation parameter. We 
expect that the understanding of the symmetry properties of prefect-fluid solutions will help 
to incorporate later the rotation parameter. 

Although there exists quite a large number of exact solutions JfJ, only a few can be 
considered as physically meaningful. This type of solutions are usually classihed in terms 
of their symmetry. Consider, for instance, spherically symmetric spacetimes. According to 
Birkhoff’s theorem, the exterior held is uniquely determined by the Schwarzschild solution. 
For the interior held, however, a quite large number of exact solutions are known. In fact, 
the explicit form of the interior solutions depends on the model used to determine the 
energy-momentum tensor of the source. In addition, if we limit ourselves to the case of 
perfect-huid sources, the search for exact solutions requires the knowledge of an additional 
equation which is usually taken as the equation of state of the fluid. The explicit form of 
the interior solution depends heavily on the properties of the equation of state. Recently, 
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several generating techniques have been proposed that allow one to obtain all spherically 
symmetric perfect-fluid interior solutions with and without equations of state 3-6]. 

In the case of stationary axisymmetric spacetimes^ only a few solutions are known for 
the interior of a rigidly rotating perfect fluid. In 18-20] several exact solutions were found 
which, however, are characterized by negative pressures. In ]21, 22] more realistic interior so¬ 
lutions were obtained with physically plausible equations of state for rigidly rotating sources. 
Nevertheless, the problem of finding the solution for the corresponding exterior held remains 
open because the attempts to solve directly the resulting vacuum held equations, together 
with the energy conditions and the matching conditions, have not been very successful. In 
view of this situation, it seems reasonable to explore other methods. In particular, we be¬ 
lieve that the methods for generating new solutions from known ones could give some new 
insight into the problematic j^]. The present work can be considered as a hrst step in this 
direction. 

In this work, we focus on non-rotating perfect-huid configurations. We propose a line 
element which is specially adapted to the investigation of the symmetry properties of the 
corresponding held equations. In fact, we will show that a particular symmetry corresponds 
to a transformation that can be used to generate new exact solutions. This paper is organized 
as follows. In Sec. El we propose for static axisymmetric interior spacetimes a particular 
metric whose held equations can be written in a relatively compact form that resembles 
the vacuum case. In Sec. m we prove a theorem about a transformation that can be 
used to generate new exact solutions, and can interpreted as a generalization of the Zipoy- 
Voorhees transformation for the vacuum limiting case. In Sec. IIV1 we present some explicit 
examples of the procedure proposed here to generate new perfect-huid solutions of Einstein’s 
equations. Finally, in Sec. IVJ we discuss our results and comment some open problems. 


II. LINE ELEMENT AND FIELD EQUATIONS 

We will consider a static axially symmetric spacetime with Killing vectors = d t and 
= dip, where t is the time-coordinate and ip is the azimuthal angle. This implies that 
dg a fj/dt = dg a p/d(p = 0. For the remaining spatial coordinates we choose 6 as the polar 
angle and r as a radial-like coordinate. As for the line element for this type of spacetimes, 
several choices are possible Q. In the case of interior solutions, a particular choice [7] has 
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been intensively used to derive and analyze approximate numerical solutions. In principle, 
for a given set of symmetry conditions and coordinates, all possible line elements should be 
equivalent. However, the point is that the structure of the held equations depends on the 
explicit form of the line element and, therefore, a particular choice might be more suitable 
for the investigation of the held equations structure. In this work, we propose the following 
line element 


ds 2 


e 2 ^dt 2 - e~ 2 ^ 




+ /. i 2 d(f 2 , 


( 1 ) 


to study the structure of the held equations with a perfect-huid source. Here ip = 

7 = q(r, 9), p = /i(r, 9), and h = h(r). Moreover, we use geometric units with G — c — 1. 
The Einstein equations for a perfect fluid with 4-velocity U a , density p, and pressure p 


R a p ~ 7 ^ R 9 ap = 87 r [(p + p)U a Up - pg af) ] 


( 2 ) 


for the line element (JT]) lead to a system of four independent differential equations. The 
main equations, relating the metric functions //, ip, and h, can be written as 
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where the comma represents partial dif ferentiation with respect to the corresponding co¬ 
ordinates, and we have introduced the “normalized” density p and pressure p by means 
of 

( 5 ) 


p = pe 2 ^ 


p = pe 2 ' 7 


Notice that there is no independent equation for the metric function h(r). The reason is 
that in the line element (ITJ) it is possible to introduce a new radial-like coordinate, say 
dR = dr/yjh , which absorbs the function h, implying that it does not represent a genuine 
degree of freedom. Nevertheless, we keep h as an auxiliary function which can be chosen 
arbitrarily, in principle. Accordingly, if we attempt to solve the above system it is necessary 
to specify h, by using some additional criteria. This will turn to be useful for the purpose 
of the present work, as we will show below. As for the metric function 7 , it is determined 
by two first order differential equations 
1 

1,r T o i 2~ 

h 9r + 9,6 


9 


9,i 






^P,6^P,r P,rP p,Qp,rO P,rP,00 


( 6 ) 


4 













1 


7,0 = / 7/i 2 + ^2 AM (V >?0 - ^ 2 r) + 2hn } r4>,04>, r ~ 8 tTP y6 P + hfl^ r g + AM/i^j , (7) 

which can be integrated by quadratures once the main held equations (J3J) and (JH) are solved, 
and the normalized pressure p is known. It is important to notice that if we introduce 
the differential equations (l3l)-(ITj) into the original Einstein equations (J2]), a second order 
differential equation for 7 emerges 


, 7,00 , / 2 , Vi 2 

7 ,rr H-^-1" r,r + 


0 , ^,r7,r 

7T 2h 


87T _ 

= T P ’ 


( 8 ) 


which must also be satisfied. Nevertheless, it can be shown that this equation is identically 
satisfied if the two first order differential equations (EJ) and (J7J) for 7 and the conservation 
equation for the non-normalized parameters of the perfect fluid 


P,r = ~{p + P)ip,r , P,e = ~(p + p)i/),o (9) 

are satisfied. 

Notice that the structure of the differential equations as presented above resembles the 
structure of the vacuum field equations. In fact, for vanishing density and pressure with 
h = 1 , the equation for /i is trivially satisfied, if /i is used as a radial-like coordinate (p ~ r), 
and the only remaining field equation is (@J for the function -0 whose solution determines 
uniquely the remaining function 7. Since the conservation law is identically satisfied in this 
case, the compatibility condition following from Eqs. () 6 |) -(l 8 |) coincides with the main field 
equation (|4]). 

Notice also that the present choice of the line element leads to a particularly simple gen¬ 
eralization of the Tolman-Oppenhcimer-Volkov (TOV) equation for spherically symmetric 
perfect-fluid sources. This has been also obtained previously with different choices of line 
elements [7]. In fact, the TOV equation coincides with the first equation given in ([9]), with 
r being the radial spherical coordinate. The generalization to include the axisymmetric case 
implies a second TOV equation for the angle coordinate 6. This simple form of the TOV 
equations allows us to easily calculate the pressure, if p and 0 are given. 


III. THE TRANSFORMATION 

The search for physically reasonable exact solutions of the above system of differential 
equations is in general a very complicated task. In principle, we have a system of four 
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independent equations ©-([7J) for the six unknowns p, 0 , 7, h, p, and p. To close the 
system it is necessary to specify h = h(r) and, for instance, an equation of state of the form 
V — p(p)- We investigated the simple case of a polynomial function h{r ) and a barotropic 
equation of state p = up , where u is the barotropic constant factor, and found several 
approximate solutions which, however, turned out to be characterized by either negative 
values of the pressure or singularities at the origin. Although this type of solutions could 
have some mathematical interest, their physical interpretation is difficult due to the presence 
of singularities, negative pressures or negative energy densities. 

However, we investigated the structure of the field equations and found certain symmetries 
that allow us to generate new solutions from a given seed solution. This generating technique 
is based upon a transformation of the functions and coordinates entering the field equations. 
The result can be formulated as follows: 

Theorem: Let an exact interior solution of Einstein’s equations ©-© for the static 
axisymmetric line element (JT]) be given explicitly by means of the functions 

h 0 = ho(r ), p 0 = Po(r, 0 ), A = A(r, 0), 70 = 7o(a 9), Po = Po(r, 9 ), p 0 = p 0 (r, 9) . (10) 
Then, for any arbitrary real value of the constant parameter 5 the new functions 

h = ho(r), p = po(r, 9), 0 = 50 o (r, 9), p = Sp 0 (r , 9), p = 5p 0 (r, 9), 9 = -^= , (11) 


7 (r, 9) = S 2 7 o(r, 9) + (<5 2 - 1) 


v~ a 


ho + v 2 


dr 


ro 


+87t5(1 — S) 


M 0 ,r^° A , h0,f) + 

- t dr + k, v = —— , k = const. 


ho + Po,r 

represent a new solution of the field equations 0-0. 

Proof: Since the set of functions (fTOl) represents a solution, the field equations 
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(14) 


ho \2/i(j po 

are identically satished. Notice that the normalized quantities po and po, as written in the 
above identities, are just explicitly known functions of their arguments so that we can treat 
them as independent functions. In other words, at the level of the identities, which follow 
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from the field equations for the seed functions, we can treat po, po, Vh and 70 as algebraically 
independent functions. 

If we now replace the functions h = ho/8 and p = fj, 0 (r,6) in Eq. (131) . and perform the 
coordinate reparametrization 9 —* 9 = 9 /\f8, we obtain the equation 

167T 


/^O ,rr ho ,99 ho, r po, t 


-5p . 


(15) 


Po hopo 2/io/io ho 
This equation coincides with the identity (ITS]) only if 9 —> 9 and 8 p p 0 . So, in general, it 
is not identically satisfied. However, one can interpret it as the definition of a new pressure 
Pnew(r,9) = 8p{r,6). To establish the connection with the seed solution we assume that 
Pnew( r ,9) = 8po{r,9). It then follows that the functions h = h 0 and p = p, 0 (r,9) satisfy the 
field equation (J3]) for a new pressure 8 po(r, 9). The next step is to show that these functions 
can be made to satisfy the remaining field equations. 

Consider the field equation (JTJ) for if)(r,9) with p = p 0 (r, 9), h = h 0 (r ) and the new 
pressure p{r,9 ) = 8po(r,9). Then, we obtain 


w-++ 

ho 


ho,r ^ PO,r 


1 , ho,e^e 47T 
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(16) 


„ 2ho Po 

A simple inspection of this equation shows that it reduces to the identity (TT4|) for the choice 
ip = 8ipo(r, 9) and p = 8p 0 (r, 9). This proves that the functions (III]) determine a new solution 
of the field equations (J3]) and (J4j). 

We now calculate the explicit form of the function 7 . First of all, we notice that since 
7 o(r, 9) is a solution so is also 70 (r, 9) when we change 9 by 9 everywhere in the equations. 
Then, the corresponding field equations become identities, i. e., 
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H(r , 0 ) = h 0 p 2 0r + Pop ■ 

We now introduce the new solution (TTTD into the field equation (j 6 ]), and obtain 

8 2 P 0 [po,r (h 0 ^ 0 , r ~ + 2 Po, 9 ^ 0 , 9 ^ 

-\~8tt 5 poPo,rPo + PopPo,r9 — PO,rPo,99 


99 
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Using Eq. fjTTj) to replace the term in squared brackets of the last equation, after some 
algebraic rearrangements, we get 

d 2 - 1 ' ‘ 


7,r ~ <5 70 ,r + 


Vo ,0 
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Consider now the second equation (JTJ) for the function 7 . Introducing the values of the 

new solution (TITj) into Eq.fJTj), we obtain 
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7.1 = 


H(r, 0) 


^Vo [Vo,g 0,0 “ h oi’lr) + 2h 0V O.r^o.^0, 

StyS fXQfiQ Qpo + (hoVo,rVo,r6 Vo,oVo ,01 


( 22 ) 


Furthermore, we use the identity (fTHjl to replace the expression contained in the squared 
brackets. It is then straightforward to get 


7 , 0 - = £ 2 7 o , 0 ~ + ^(1 - 52 )^ E lnH (v 0 ) + 8 ^ 5(5 - i )^ 0 -^ 0 


(23) 


2 V " ' 7 90 v, ' / ' ""' v ' ' H( rj 6) 

The above equation can immediately be integrated, and the solution is determined up to an 

arbitrary function of r 

1 (r,e) = S 2 l0 + kl-S 2 )lnH + SirS{S-l) f^^.S+F(r) (24) 

which can be fixed by considering the second equation (l2Tj) for the function y(r, 6). We 
obtain 

1 

, ^^dr - - 

hoVo,r + Vo,e 2 
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where k is an integration constant. Inserting this expression into Eq. (j24lh after some alge¬ 
braic rearrangements we finally obtain 

tfn*)-*Vn*) + (*-i), fto + [/2 - ■ . J um 

which is equivalent to Eq. (TT 2 lh This ends the proof of the theorem. 
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U ° dr + 87 t 5(1 — S) f ^°^ 0 ’ r ?° dr + k , (26) 


IV. EXAMPLES 

The theorem proved in the previous section states that given a static axisymmetric so¬ 
lution it is possible to generate a new solution which contains an additional new parameter 
5. In this section, we will apply the transformation to show that it can be used to generate 
new solutions that are physically different from the seed solution. 



















A. The vacuum q-metric 


Let us consider the vacuum limiting case of the above transformation. The particular 
choice 

p — 0, p = 0, /i = r, h — 1, d = z (27) 

in Eq. CQ) leads to the static axisymmetric Lewis-Papapetrou line element for vacuum space- 
times 

ds 2 = e^dt 2 - e~ 2 ^ e 2l (dr 2 + dz 2 ) + r 2 dip 2 (28) 

in cylindrical coordinates ( t,r,z,tp ) for which the held equations (J3])-([7|) reduce to 


Vv + -f/v + ”0,22 = 0 , 7,r = r(0 2 - 0 2 ), 7* = 2? y/y02 , 


(29) 


for the functions 0 and 7 only. Then, according to the Theorem proved in the previous 
section, if 0 O = tM'L z) and 70 = 70 (r, - 2 ) represent a particular solution of Eqs. fl29l) . a new 
solution 0 , 7 is given by 


= Sij> 0 (r t z), 7 (r,z) = * 5 2 7 0 (r,z) . 


(30) 


Notice that for this new solution the line element corresponds to with z —> z. Since 
the new functions depend also on z, one can use the coordinate 0 instead, without loss o 
generality. The transformation ()30l) was originally proposed by Zipoy js] and Voorhees 
and has been investigated in several works (see, for instance, 10-|17j and references therein). 
For this reason, the transformation defined in Eqs. fllip and m can be considered as a 
generalization of the Zipoy-Voorhees transformation that includes the case of perfect-fluid 
spacetimes. 

A particular example of the Zipoy-Voorhees transformation arises if we take the 
Schwarzschild metric as seed solution as follows. First, we introduce spherical coordinates 
(f, r, d, p) into the line element (T28|) by means of the transformations 


r = m + -(r + + r -) > cost? = (r + — r_) , r± = \Jr 2 + (z ± m) 2 


(31) 


and obtain 


ds 2 = e^dt 2 - e~W 


( 2m m 2 sin 2 tP 

e 2 M 1 - — + --- 


dr 2 




2 _ 2m 


+ r 2 dd 2 


+ f 2 (l-^)sm 2 *V 


(32) 
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Then, the Schwarzschild metric corresponds to the particular solution 


1 r — 2 m 

^0 = 7 ln -~— 

2 r 


1 r 2 — 2 mr 

2 f 2 — 2 mr + m 2 sin 2 § 


(33) 


Applying the transformation (1301) with 5 = 1 + q, i.e., 


, 1 ,, , . f — 2m 
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the resulting line element can be expressed as 
2 m \ 1+9 


ds = 


r / 

I 2m \ -9 

f / 


dt~ 


1 + 


m 2 sin 2 i?\ 9(2+9) / dr 2 


r 2 — 2 mr 


| _ 2 m 


+ r 2 cm 2 + r 2 sin 2 '/My ? 2 


,(35) 


which, as expected, reduces to the Schwarzschild metric in the limiting case q —> 0. 


A detailed analysis of this metric shows that m and q are constant parameters that 
determine the total mass and the quadrupole moment of the gravitational source [rij. The 
metric (1351) has been interpreted as the simplest generalization of the Schwarzschild metric 
with a quadrupole. I11 the literature, this metric is known as the 6— metric, referring to 
the parameter <5, or as the 7 —metric for notational reasons. We use the term q —metric 
to emphasize the physical interpretation of the new solution in terms of the quadrupole 
moment. 

Whereas the seed metric is the spherically symmetric Schwarzschild solution, which de¬ 
scribes the gravitational field of a black hole, the generated q —metric is axially symmetric 
and describes the exterior field of a naked singularity 17]. This shows that the transfor¬ 
mation presented in the previous section generates non-trivial solutions even in the case of 
vacuum solutions. 


B. Generalization of the interior Schwarzschild solution 


One of the most important interior solutions of Einstein’s equations is the spherically 
symmetric Schwarzschild solution which describes the interior field of a perfect-fluid sphere 
of radius R and total mass m. The corresponding line element can be written as 


ds~ = 


:f(R) ~ o/( r ) 


dt 2 


dr 2 

f 2 ( r ) 


r 2 (d6 2 + sin 2 9d(p 2 ) , 


(36) 
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with 


f(r) = V 1 “ 


2 mr 2 
i ? 3 


(37) 


The parameters of the perfect fluid are the constant density po and the pressure po which is 
a function of the radial coordinate r only 

m - m 


Po = Po: 


(38) 


3 f(R)-f(r) 

We will use the above metric as the seed solution (fTUD for the general transformation 
(fill) . Then, a straightforward comparison with the general line element (JTJ) yields 

e^° = — i/(r) , ho = r 2 / 2 (r) , p 0 = rsinOe , e 7(l = re^° . (39) 

According to the Theorem presented in Sec JIIIl the new solution can be obtained from 
Eq. (139|) by multiplying the corresponding metric functions with the new parameter 5. Then, 
the new line element can be represented as 

dr 2 


ds 2 = e 2H ’°dt 2 - e~ 25 ^ 0 


=27 


r 2 y 2 ^ 


+ dd 2 + r 2 e 2 ^° sin 2 ddp 2 


(40) 


where the new function 7 is given by 

(1 — 5 2 ) + 87 t 5(1 — 5) sin 2 9r 2 p 0 


7 = ^ 7o + 


r/ 2 (r)(l + ripor) sin^ 9 + 


l+ripO‘. 


cos 2 6 


dr + k 


where 


i’or = 


2 mr 


i? 3 /(r)[3/(i?)-/(r)] ‘ 
Moreover, the parameters of the perfect fluid source are 


(41) 


(42) 


p = 3p 0 e 2 T«- 2 '>+ 2 ( J - 1 W’» ,p = 6p„e' 2 ' , °- 2 ~< +2(s - 1) *° . 


(43) 


Notice that the new function 7 as given in Eq. (14 1 j) depends explicitly on the coordinate 9, in 
contrast to the seed metric function y 0 which depends on the coordinate r only. This proves 
that the new solution is not spherically symmetric, indicating that it is physically different 
from the interior Schwarzschild solution. This conclusion is corroborated also by the fact 
that the density and pressure of the new solution are functions of the angular coordinate 
too. 
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V. CONCLUSIONS 


In this work, we presented a new method for generating perfect-fluid solutions of the 
Einstein equations, starting from a given seed solution. The method is based upon the 
introduction of a new parameter at the level of the metric functions of the seed solution in 
such a way that the generated new solution is characterized by physical properties which are 
different from those of the seed solutions. This means that the solution generating technique 
proposed in this work can be used to generate non-trivial new solutions. 

To show the validity of the method we propose a line element which is especially adapted 
to handle the problem. In fact, the perfect-fluid field equations turn out to be split into a set 
of two partial differential equations that once solved can be used to integrate by quadratures 
the remaining equations. We also introduce in the line element an auxiliary function that 
is not fixed by the field equations and, therefore, can be used to simplify the analysis of the 
main field equations. 

The method has a well-defined vacuum limit. In fact, for a particular choice of the 
metric functions and for vanishing density and pressure, we obtain the Lewis-Papapetrou 
line element for static axisymmetric vacuum fields. In this limit, the generating technique 
proposed here coincides with the Zipov-Voorhees transformation which can be used to gen¬ 
erate vacuum solutions from vacuum solutions. As a particular example we derive from the 
Schwarzschild metric a solution with quadrupole moment, which we call the q— metric, and 
can be interpreted as describing the simplest generalization of the Schwarzschild spacetime 
with a quadrupole parameter. 

In the case of perfect-fluid solutions, we use the interior spherically symmetric 
Schwarzschild solution with constant density to generate a new interior solution which turns 
out to be axially symmetric. 

The generation of further new solutions does not represent any particular difficulty for 
the method developed in this work. The important problem of matching the so generated 
interior solutions with the corresponding exterior solutions is beyond the scope of the present 
work, and will be the subject of future investigations. 
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